A  UNIFIED  MODEL  FOR  HYSTERESIS 
IN  FERROIC  MATERIALS 


Ralph  C.  Smith 

Department  of  Mathematics 
North  Carolina  State  University 

•  Marcelo  Dapino,  Ohio  State  University 

•  Zoubeida  Ounaies,  Virginia  Commonwealth  University 

•  Stefan  Seelecke,  North  Carolina  State  University 


Research  Support: 

Air  Force  grant  AFOSR-F49620-01- 1-0107 
NASA  grant  NAG-1-01041 
NSF  grant  CMS-0099764 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

00  JUN  2003 

2.  REPORT  TYPE 

N/A 

3.  DATES  COVERED 

4.  TITLE  AND  SUBTITLE 

A  Unified  Model  For  Hysteresis  In  Ferroic  Materials 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROIECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Department  of  Mathematics  North  Carolina  State  University;  Ohio  State 
University 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

See  also  ADM001697,  ARO-44924.1-EG-CF,  International  Conference  on  Intelligent  Materials  (5th) 

(Smart  Systems  &  Nanotechnology).,  The  original  document  contains  color  images. 

14.  ABSTRACT 

15.  SUBIECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

uu 

18.  NUMBER 
OF  PAGES 

17 

19a.  NAME  OF 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


TALK  OVERVIEW 


I.  Smart  Material  Applications 

•  Piezoelectrics  -  Tunable  lenses,  nanopositioning 

•  Shape  Memory  Alloys  -  Shape  modification 

•  Magnetostrictives  -  High  force  transduction 

II.  Modeling  Hierarchies 

•  Develop  energy  relations  at  lattice  level  for  single  crystal, 
homogeneous  compounds. 

•  Incorporate  polycrystallinity  and  material  nonhomogeneities 
through  stochastic  homogenization  to  obtain  macroscopic 
constitutive  relations. 

III.  Experimental  Validation 

•  PZT5A  -  Major  and  biased  minor  loops 

•  Terfenol-D  -  Major  and  biased  minor  loops 

•  SMA  -  Thin  film 

Goal: 

•  Develop  a  unified  framework  for  model  and  control 
design  for  ferroic  (ferromagnetic,  ferroelectric  and 
ferroelastic)  materials. 


I.  PIEZOCERAMICS  -  NANOPOSITIONING 


Atomic  Force  Microscope: 

•  Can  provide  atomic  resolution 

•  Human  chromosomes 
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Model  Development: 

•  Must  be  able  to  attain 
Angstrom-level  resolution 

•  Must  be  implemented 
in  real  time 

Future  Capabilities: 


Electric  Field  (V/m) 


•  Nanopositioning  at  the  atomic  level 


•  Electron  spin  detection:  requires  angstrom-level  tracking 


SHAPE  MEMORY  ALLOYS 


Aerospace  and  Aeronautic  Applications: 

•  Shape  modification 
of  airfoils 

•  Construction  of  modular 
mirrors  and  antennas 

•  Future  Directions:  MEM’s, 
thin  films  -  Sandia  National  Labs 


Embedded  SMA 
Actuators 


Vibration  Attenuation: 


Membrane 

Mirror 


•  Vibration  attenuation  in  aerospace  structures  (e.g.,  mirrors) 

•  Operate  in  hysteretic  regime  to  optimize  damping 


MAGNETOSTRICTIVE  MATERIALS 


Automotive  Application:  e  g.,  High  speed  milling 


Aerospace  Applications: 

•  Magnetostrictive  composites  for  active  structural  damping 
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MATERIAL  ANALOGIES 


Ferroelectric 

Ferromagnetic 

Ferroelastic 

Polarization 

Magnetization 

Strain 

Electric  field 

Magnetic  field 

Stress 

Paraelectric  phase 

Paramagnetic  phase 

Austenite  phase 

Ferroelectric  phase 

Ferromagnetic  phase 

Martensite  phase 

Ferroelectric 
domain  walls 

Bloch  or  Neel  walls 

Boundaries  between 
martensite  variants 

Devonshire  theory 

Mean  field  (Weiss)  theory 

Micromagnetic  theory 

Landau  theory 

Ginzburg-Landau  theory 
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MECHANISMS  WHICH  PRODUCE  HYSTERESIS 

Ferroelectric:  Dipole  switching 


Ferromagnetic:  Moment  rotation,  domain  wall  losses 
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MESOSCOPIC  MODEL  -  HELMHOLTZ  ENERGY 


Helmholtz  Energy: 

i/>(P,T)  =  U  -  ST 


$oiV[l  -  ( P/Ps )2]  +  TkN 


4V 
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Plnl^P)+Pln(l-(P/P'2 


Helmholtz  Energy  Relations: 
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MESOSCOPIC  MODEL  -  GIBBS  ENERGY 


V(P,T)  =  G(0,P,T) 
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Gibbs  Energy  Relations: 

G{E,P,T)=^(P,T)-EP 

G(H ,  M,  T)  =  ^(M,  T)  - 
G(cr,  £,  T)  =  T)  —  as 


(Ferroelectric) 

(Ferromagnetic) 

(Ferroelastic) 


MESOSCOPIC  MODELS  -  UNIFORM  LATTICE 


Boltzmann  Probability: 

/j,(G)  =  Ce-GV/kT  =>  p+ _ 


kT 


\  2irmV 2/3 


e-G(E,P0(T),T)V/kT 

JPq 


Evolution  Relations: 


dx+ 
dt 
dx _ 
dt 
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Ferroelastic:  Analogous 


Polarization: 

P  =  x+  (P+)  +  x_  (P_) 

where 

(P+)  =  /“  PM(G)dP 


Magnetization: 

M  =  x+  (M+)  +  x_  (M_) 

where 

<M+>  =  Q  Mp(G)dM 


MACROSCOPIC  MODELS 


f(Ec)  =  Cie~ME c/^c)/26]2 

Variable  Effective  Fields: 

f{£)  =  C2e-{£~E W 


Macroscopic  Models: 

[P(E,  T)](t)  =  Cf0°°  IZ  f(Ec)f(£)[P(£  +  E,  T)](t)d£dEc 
[M(H,  T)](t)  =  C  /0°° IZ  f(Hc)f(H)[M(n  +  H,  T)](t)dHdHc: 


=  C  /0°°  IZ  f(S)f{aE)[e{aE  +  a,T)](t)daEdS 


III.  EXPERIMENTAL  VALIDATION 


PZT5A:  Lognormal/normal  densities 


PZT5A:  ID  of  general  densities  [B.  Mukherjee,  S-F.  Liu  -  Data] 


VALIDATION  -  FERROMAGNETIC  MATERIALS 

Example  1:  Steel  data  from  [Jiles  and  Atherton;  1984] 


Example  2:  Terfenol-D  Transducer 


EXPERIMENTAL  VALIDATION  -  SMA 


Thin  Films: 


CONCLUDING  REMARKS 


Model  Attributes: 

(1)  Technique  provides  unified 
framework  for  modeling  hys¬ 
teresis  in  ferroelectric,  ferro¬ 
magnetic  and  ferroelastic 
materials. 

(2)  Low  number  (5-10)  and 
physical  nature  of  para¬ 
meters  facilitates  parameter 
estimation  and  control  implemen 
mentation. 

(3)  Model  guarantees  closure 
of  biased  minor  loops. 


Field 


(4)  Model  is  amenable  to  inver¬ 
sion  which  facilitates  linear 
control  design. 

(5)  Method  provides  energy 
basis  for  Preisach  models. 
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